********************************************************
***5 January 2022
***replication files for
***Timothy Hellwig & Dani M. Marinova, "Evaluating the Unequal Economy: Poverty Risk, Economic Indicators, and the Perception Gap"
***forthcoming, POLITICAL RESEARCH QUARTERLY
********************************************************
***the following first calls up data file "Hellwig_Marinova_PRQ_EES.dta" to produce Tables 1 and 2 and Figure 1 in manuscript and appendix Tables A5-A16 and Figure A2
***it then calls up data file "Hellwig_Marinova_PRQ_EB.dta" to produce Table 3 in the manuscript and appendix Tables A17-A2 
***send questions or report bugs to thellwig@indiana.edu
********************************************************

use "Hellwig_Marinova_PRQ_EES.dta", clear

/*TABLE 1*/
meologit nret1234 pov_riskt1 i.year [pweight=1/stdp] || pov_groups : 
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_pov00 i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :

/*TABLE 2*/
meologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
meologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
meologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year [pweight=1/stdp] || pov_groups :, 

/*FIGURE 1*/
qui meologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, at(chng_gdptt1=2.40 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_une=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("A. Growth Increase") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(growth_incr, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, at(chng_gdptt1=-7.66 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_une=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("D. Growth Decrease") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(growth_decr, replace)

qui meologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, at(chng_une=-1.63 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("B. Unemployment Decrease") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(une_decr, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, at(chng_une=3.373 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("E. Unemployment Increase") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(une_incr, replace)

qui meologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year [pweight=1/stdp] || pov_groups :, 
margins, at(chng_pov00=-4.19 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_une=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("C. Poverty Rate Decrease") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(pov_decr, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year [pweight=1/stdp] || pov_groups :, 
margins, at(chng_pov00=3.21 pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_une=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("F. Poverty Rate Increase") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(pov_incr, replace)

graph combine growth_incr une_decr pov_decr growth_decr une_incr pov_incr, scheme(s1mono) ycommon
graph save Graph "figure1.gph", replace


/*********************************/
/***********APPENDIX**************/
/*********************************/
/*Tables A1-A4 & Figure A1 are produced using EU-SILC micro-level cross-section data*/
/*Replication data not available due to regulations governing use of SILC data (https://ec.europa.eu/eurostat/web/microdata/european-union-statistics-on-income-and-living-conditions)*/

/*TABLE A5*/
meologit nret1234 pov_risk i.year || pov_groups : 
meologit nret1234 pov_risk our0813 voted_pm university unemployed i.class_selfplace i.news i.year || pov_groups : 
meologit nret1234 pov_risk our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 i.year || pov_groups : 
meologit nret1234 pov_risk our0813 voted_pm university unemployed i.class_selfplace i.news chng_une i.year || pov_groups : 
meologit nret1234 pov_risk our0813 voted_pm university unemployed i.class_selfplace i.news chng_pov00 i.year || pov_groups : 
meologit nret1234 pov_risk our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une chng_pov00 i.year || pov_groups : 

/*TABLE A6*/
meologit nret1234 c.pov_risk##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year || pov_groups : 
meologit nret1234 c.pov_risk##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year || pov_groups : 
meologit nret1234 c.pov_risk##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year || pov_groups : 

/*TABLE A7 MODELS 1 & 2*/
gen subjliv = q120
recode subjliv 77/max=.
label variable subjliv "subj standard living 2009 only"
mixed subjliv pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news i.year || pov_groups : 
meologit diff_bills pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news || pov_groups : 

/*TABLE A8*/
gen subjinc1 = (subjliv-1)/6
recode subjinc1 .=0
gen subjinc2 = (social_standing-1)/9
recode subjinc2 .=0
gen subjinc = subjinc1 + subjinc2
recode subjinc *=. if q120==77 | q120==88 | d61==-8
gen subjinc2009 = (subjliv-1)/6
pwcorr pov_riskt1 our0813 ca scp msf lsf bc othrocc class_selfplace subjinc2009 university unemployed, obs sig

/*TABLE A9*/
qui meologit nret1234 pov_riskt1 ca scp bc lsf our0813 voted_pm university unemployed i.class_selfplace i.news subjinc2009 [pweight=1/stdp] || pov_groups :
predict corr if e(sample)
meologit nret1234 ca scp bc lsf our0813 voted_pm university unemployed i.class_selfplace i.news subjinc2009 [pweight=1/stdp] if corr~=. || pov_groups :
meologit nret1234 pov_riskt1 i.year [pweight=1/stdp] if corr~=. || pov_groups :
meologit nret1234 pov_riskt1 ca scp bc lsf our0813 voted_pm university unemployed i.class_selfplace i.news subjinc2009 [pweight=1/stdp] if corr~=. || pov_groups :

/*TABLE A10*/ 
meologit nret1234 pov_riskt1 pol_int female age_u40 othrnatl ca scp bc lsf othrocc i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 pol_int female age_u40 othrnatl ca scp bc lsf othrocc our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 pol_int female age_u40 othrnatl ca scp bc lsf othrocc our0813 voted_pm university unemployed i.class_selfplace i.news chng_une i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 pol_int female age_u40 othrnatl ca scp bc lsf othrocc our0813 voted_pm university unemployed i.class_selfplace i.news chng_pov00 i.year [pweight=1/stdp] || pov_groups :
meologit nret1234 pov_riskt1 pol_int female age_u40 othrnatl ca scp bc lsf othrocc our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :

/*TABLE A11*/
ologit nret1234 i.ccode pov_riskt1 i.year [pweight=1/stdp]
ologit nret1234 i.ccode pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news i.year [pweight=1/stdp] 
ologit nret1234 i.ccode pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 i.year [pweight=1/stdp] 
ologit nret1234 i.ccode pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une i.year [pweight=1/stdp] 
ologit nret1234 i.ccode pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_pov00 i.year [pweight=1/stdp]
ologit nret1234 i.ccode pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une chng_pov00 i.year [pweight=1/stdp] 

/*TABLE A12*/
ologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp]
ologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] 
ologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year [pweight=1/stdp]  

/*TABLE A13*/
mixed nret_pos pov_riskt1 i.year [pweight=1/stdp] || pov_groups : 
mixed nret_pos pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news i.year [pweight=1/stdp] || pov_groups :
mixed nret_pos pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 i.year [pweight=1/stdp] || pov_groups :
mixed nret_pos pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une i.year [pweight=1/stdp] || pov_groups :
mixed nret_pos pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_pov00 i.year [pweight=1/stdp] || pov_groups :
mixed nret_pos pov_riskt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :

/*TABLE A14*/
mixed nret_pos c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
mixed nret_pos c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
mixed nret_pos c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news c.chng_gdptt1 chng_une i.year [pweight=1/stdp] || pov_groups :, 

/*TABLE A15*/
meologit nret1234 pov_riskt1 i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups : 
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.news i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.news chng_gdptt1 i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.news chng_une i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.news chng_pov00 i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups :
meologit nret1234 pov_riskt1 our0813 voted_pm university unemployed i.news chng_gdptt1 chng_une chng_pov00 i.year [pweight=1/stdp] if class_selfplace>1 || pov_groups :

/*TABLE A16*/
/***note: not able to estimate with pweights*/ 
***1. Test the proportional odds assumption with omodel logit***
***A negative value indicates that the effects are not equal across categories (this is the case)***
tab year, gen(year_)
tab class_selfplace, gen(class_selfplace_)
tab news, gen(news_)
omodel logit nret1234 pov_riskt1 year_2 
omodel logit nret1234 pov_riskt1 our0813 voted_pm university unemployed class_selfplace_2 class_selfplace_3 news_2 year_2  
omodel logit nret1234 pov_riskt1 our0813 voted_pm university unemployed class_selfplace_2 class_selfplace_3 news_2 chng_gdptt1 chng_une chng_pov00 year_2 
***2. Fit a generalized linear model that relaxes the proportional odds assumption***
gologit2 nret1234 pov_riskt1 year_2
gologit2 nret1234 pov_riskt1 our0813 voted_pm university unemployed class_selfplace_2 class_selfplace_3 news_2 year_2  
gologit2 nret1234 pov_riskt1 our0813 voted_pm university unemployed class_selfplace_2 class_selfplace_3 news_2 chng_gdptt1 chng_une chng_pov00 year_2 

/*FIGURE A2*/
qui meologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, dydx(chng_gdptt1) at(pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_une=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("A. ME of Growth") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_growth, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_gdptt1 our0813 voted_pm university unemployed i.class_selfplace i.news chng_une chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, dydx(pov_riskt1) at(chng_gdptt1=(-10(2.5)5) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_une=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("D. ME of Poverty risk") xtitle("Growth") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_prg, replace)

qui meologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, dydx(chng_une) at(pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("B. ME of Unemployment") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_une, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_une our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_pov00 i.year [pweight=1/stdp] || pov_groups :,
margins, dydx(pov_riskt1) at(chng_une=(-2(2)6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_pov00=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("E. ME of Poverty risk") xtitle("Unemployment") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_pru, replace)
	
qui meologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une i.year [pweight=1/stdp] || pov_groups :,
margins, dydx(chng_pov00) at(pov_riskt1=(0(.1).6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_une=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("C. ME of Poverty Rate") xtitle("Poverty risk") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_pov, replace)
qui meologit nret1234 c.pov_riskt1##c.chng_pov00 our0813 voted_pm university unemployed i.class_selfplace i.news chng_gdptt1 chng_une i.year || b: pov_riskt1 || pov_groups :, 
margins, dydx(pov_riskt1) at(chng_pov00=(-6(2)6) our0813=.08 voted_pm=0 university=0 unemployed=0 class_selfplace=2 news=1 chng_gdptt1=0 chng_une=0 year=2014) pr(outcome(4)) pr(outcome(1)) post vsquish
marginsplot, recast(line) recastci(rarea) graphregion(fcolor(white) ifcolor(white)) ///
	title("F. ME of Poverty risk") xtitle("Poverty rate") ytitle("Econonomic evaluation") ///
	yline(0) scheme(s1mono)	ciopt(color(%20)) name(me_prp, replace)

graph combine me_growth me_une me_pov me_prg me_pru me_prp, scheme(s1mono) r(2) ycommon
graph save Graph "figure2.gph", replace

clear

use "Hellwig_Marinova_PRQ_EB.dta", clear

/*TABLE 3*/
melogit gdp_correct1 c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
melogit une_correct1 c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
melogit povabs_correct c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
	
/*TABLE A7 MODEL 3*/
meologit endsmeet pov_risk0814 i.education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:

/*TABLE A17*/
pwcorr gdp_correct1 une_correct1 povabs_correct pov_risk0814 education2 une_risk0814 hh_children unemployed, obs sig

/*TABLE A18*/
melogit gdp_correct1 c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed wealth [pweight=1/stdp] || pov_groups:
melogit une_correct1 c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed wealth [pweight=1/stdp] || pov_groups:
melogit povabs_correct c.pov_risk0814 education2 une_risk0814 i.hh_children unemployed wealth [pweight=1/stdp] || pov_groups:
	
/*TABLE A19*/
melogit gdp_correct1 c.pov_risk0814##c.education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
melogit une_correct1 c.pov_risk0814##c.education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
melogit povabs_correct c.pov_risk0814##c.education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
	
/*TABLE A20*/
melogit gdp_correct5 pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
melogit une_correct5 pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:
gen lngdp_diff = ln(gdp_diff+1)
gen lnune_diff = ln(une_diff+1)
mixed lngdp_diff pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:, var
mixed lnune_diff pov_risk0814 education2 une_risk0814 i.hh_children unemployed [pweight=1/stdp] || pov_groups:, var


